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ABSTRACT 

Recent theoretical work has solidified the viability of the collisional runaway scenario in young dense 
star clusters for the formation of very massive stars (VMSs), which may be precursors to intermediate- 
mass black holes (IMBHs). We present first results from a numerical study of the collisional runaway 
process in dense star clusters containing primordial binaries. Stellar collisions during binary scattering 
encounters offer an alternate channel for runaway growth, somewhat independent of direct collisions 
between single stars. We find that clusters with binary fractions > 10% yield two VMSs via collisional 
runaways, presenting the exotic possibility of forming IMBH-IMBH binaries in star clusters. We 
discuss the implications for gravitational wave observations, and the impact on cluster structure. 
Subject headings: stellar dynamics — globular clusters: general — intermediate-mass black holes — 
methods: n-body simulations 


1. INTRODUCTION 

Observations hinting at the possible existence of 
intermediate-mass black holes (IMBHs) have mounted 
in recent years. Ultra-luminous X-ray sources—point 
X-ray sources with inferred luminosities > 10^® ergs“^— 
may be explained by sub-Eddington accretion onto BHs 
more massive than the maximum mass of ~ IOMq ex¬ 
pected via core collapse in main sequen ce stars, although 
viabl e alternative explanations exist llMiller V Colberil 
120041 . Similarly, the cuspy core velocity dispersion 
profiles of the globular clusters M15 and G1 may also 
be expl ained by the dynamical influence of a central 
IMBH ll van der Marel et HI l2002t iGerssen et al.1 l2002t 
iGebhardt et al.ll200,4l . although theoretical work sug¬ 
gests that the observations of M15 may be equally-well 
explained by a collection of compact stellar rem nants in 
the cores of the clusters llBaumgardt et al.ll2n0.‘lll . 

At least three distinct IMBH formation mechanisms 
have been discussed in the literature. The first, and 
possi bly simplest, is core collapse of a massive Pop HI 
star llMadau fc Reesll200111 . The very low metallicity of 
Pop HI stars {Z < 10~^Zq) allows much larger main- 
sequence stars to form, limits mass loss during stellar 
evolution, and increases the fraction of mass retained 
in the final BH (for stars more massive th an ~ 250Mq) 
iFrver et alJl200H IMiller fc Golberti 120041 . The second 
is the successive merging of stellar mass BHs via dy¬ 
namical interactions, which may occur in star clusters 
that do not reach deep core collapse before ~3Myr, 
when the most massive cluster stars have become 

BHs llKu lkarni et aLllQQJ USigurdsson fc HernauistllQQ.It 

iPort/gie^^wa^^^^^Hi^nl 1200(1: IMill'er'fcHainiit^ 

120021 lOTea^^^im2006lh The process is relatively 
inefficient—in terms of the amount of mass added to 
the growing BH per BH ejected from the cluster— 
requiring BH seeds > SOOMq to create lO^M© IMBHs 
lIGiiltekin et alJ 120041 . even when aided by the Kozai 
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mechanism llMiller fc HamiltonI l2002fl o r gravitational 
wave losses during close approaches lIGiiltekin et all 
l2?io^ . Introducing a mass spectrum for the BHs de- 
creas es the required seed mass, although growth is still 
rare llO’Learv et al ] |200a . The third is the runaway 
merging of main-sequence stars via direct physical colli¬ 
sions to form a very massive star (VMS), which may then 
collapse to form an IMBH llPortegies Zwart et alJj l99f ; 

El^i^jiijzi^ki /( al.|2 Q011IPortegies ZwartTTMcMiiiaril^OfU ; 

Giirkan et alJl2f)0^ . Recent work shows that runaway 

growth of a VMS occurs generically in clus ters with deep 
core c ollapse times shorter than ~ 3 Myr llFreitag et alJ 
1200 ^ . 

With the exception of one simulation 
llPortegies Zwart et al.l[2f)f)^ . all simulations of runaway 
collisional growth in clusters have ignored the effects 
of primordial binaries, which are known t o exist in 
cluste rs in dynamically significant numbers iisnr et alJ 
[1993). Indeed, some numerical results suggest that the 
primordial binary fraction (/{,) may have to be nearly 
100% to explain the currently observ ed binary fractions 
in cluster cores lllvanova et alJl2C)iO^ . Primordial bina¬ 
ries are an important piece of the runaway collisional 
growth puzzle, since they introduce two effects which 
may strongly affect the process. On the one hand, 
binaries generate energy via dynamical scattering 
interactions in cluster cores, supporting the core against 
deep collapse and limiting the maximum stellar density 
attai nable, and hence limi ting the direct stellar collision 
rate llHeggie fc Hutll2C)fi^ . On the other hand, stellar 
collisions are much more likely in dynamical interactions 
of binaries, since the interactions are t ypically resonant 
llBacon et al.lll994lFregeau et al.ll2004l . Since these two 
effects of primordial binaries act in opposite senses with 
respect to the collision rate, it is not clear a priori how 
they affect the collisional runaway scenario. 

Before appealing to numerical methods, however, one 
can gain insight into the effects of pri mordial binaries 
by considering the coagulatio n equation llLedl 19931 120001: 
iMalvshkin fc Go odmanll2001l) — a simp lification of which 
is presented in iFreitag et al.l ll2005bli — which describes 
the evolution of a spectrum of masses due to mergers. 
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For growth to occur in a runaway fashion, the coagula¬ 
tion equation requires that the cross section for collisions 
with the runaway object scales sufficiently rapidly with 
its mass: S'coii oc M’*, with 7 ; > 1. For single-single 
star collisions in star clusters in which the central ve¬ 
locity dispersion is less than the escape speed from the 
surface of a typical star (so that the cross section is domi¬ 
nated by gravitational focusing), this corresponds to the 
constraint R oc with a > 0 on the main-sequence 
mass-radius relationship, which is satisfied by main- 
sequence stars of any mass or metallicity. With s ome ap¬ 
proximations, the coagulat i on equation analysis gggid 
1^ iHut fc Bahcalll ll^ iSigurdsson fc Phinnevllim 

Giiltekin et alJl 200 f)ll can be applied to collisions occur- 
ring in binary scatter ing interactions. From Fig. 4 of 
IGiiltekin et alJ 11200511 . the cross section for close ap¬ 
proach distances of rmin in binary-single scattering en¬ 
counters scales as (<S'coii/7ra^)(uoo/'Cc)^ oc (?'min/a)^, with 
0-3 ^ 7 ^ 1 , where a is the binary semimajor axis, Voo is 
the relative velocity between the binary and single star 
at infinity, and Vr is the critical velocity (see eq. [ 12 ] of 
IGiiltekin et alJ[2?)0^ . Using the radius of the runaway 
star for rmin, R oc with 0.5 < /3 < 1 for the scaling 
of the mass-radius relation, and assuming that the bind¬ 
ing energy of the binary is roughly preserved during the 
encounter so that a oc M, we find: 

5eon oc m2+o(/3-i) , ( 1 ) 

with 0.3 ^ 7 1 and 0.5 ^ /3 ^ 1. The minimum of the 

exponent is ssl.5. Thus according to the coagulation 
equation, collisions induced in binary-single scattering 
interactions should yield runaway growth of a VMS. This 
result says nothing about the rate of growth of the run¬ 
away object. In other words, it is still not clear whether 
binary interactions will limit the cluster core density such 
that the runaway timescale is longer than the massive 
star main-sequence lifetime of « 3 Myr, in which case the 
process would be halted. 

Assuming that the cluster core density reached is high 
enough for the runaway to proceed, it would appear 
that a single binary is sufficient for a binary interaction- 
induced runaway to occur. This is, of course, not the 
case, since binary scattering interactions tend to destroy 
binaries. Thus we expect that for sufficiently low /&, 
the runaway will be primarily mediated by single-single 
collisions. For sufficiently large ft, the runaway will be 
primarily mediated by binary-binary interactions (the 
analysis above assumes binary-single interactions—for 
binary-binary the value of 7 will be smaller, still allowing 
a runaway by eq. ^). For intermediate fb it is possible 
that a binary interaction induced-runaway could proceed 
until the core binary population is sufficiently depleted 
that the cluster’s core collapses. If the first runaway is 
far enough from the center of the cluster when the core 
collapses, it is possible that a second runaway will be 
formed during core collapse, mediated by single-single 
collisions. The exotic possibility of forming two VMSs 
in a cluster, and thus two IMBHs, is a tantalizing one, 
with implications for gravitational wave observations and 
cluster dynamics. 

In this letter we present first results from a study of the 
runaway collisional scenario for the formation of VMSs in 
young dense clusters with primordial binaries. We briefly 
describe our numerical method, present results showing 
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Fig. 1.— Evolution of the cluster Lagrange radii (lower panel) 
in units of the initial cluster half-mass radius (left axis) and in 
parsecs (right axis), and mass of the (single) runaway for a cluster 
with 10® objects and /(, = 0.05. Time is shown both in Myr and 
in units of the central relaxation time. 


the growth of two runaways, and discuss the implications. 

2. NUMERICAL ANALYSIS 

We use our Monte Garlo cluster code to simulate the 
evol ution of young star clusters with pr i mordial bina¬ 

ries Joshi et ffi]l200nl: iFreg-eau et al.ll20(^ iGiirkan et alJ 
12004 ). This code uses the Henon method for two-body 
relaxation, and incorporates the Fewbody V-body inte¬ 
grator to perform dynamical scattering encounters of bi¬ 
naries. Stellar collisions are handled by assuming that 
stars whose surfaces touch merge with no mass loss, an 
assumption that has been shown to be valid for clusters 
with low veloci ty dispersion, such as globulars or young 
dense clusters llFreitaer et al.ll20n5blar . Gollisions are al¬ 
lowed to occur directly in singl e-single star encoun ters, 
and during binary interactions llFregeau et alJl2?)03) . 

For initial conditions we assume a Wq = 3 King model 
density profile for the cluster, with no initial mass seg¬ 
regation. All stars are main-sequence stars with masses 
in the range 0.2 < M/Mq < 120, distributed accord¬ 
ing to a Salpeter mass function. Simulations show that 
the primary condition required for a runaway is that the 
core collapse time is shorter than the main -sequence life¬ 
time f or the most massive stars, 3 Myr llFreitag et alJ 
1200 ^ : and that for clusters of single stars with a 
wide mass spectrum, the core collapse time is always 
tec ~ 0.15trc(0), (where trc(O) is the initial relaxation 
time in the cor e), independent of c luster mass, size, or 
density profile (IGiirkan et alJ 1^041) . This allows us to 
set the central density of the cluster such that the pre¬ 
dicted core collapse time is either less or greater than 
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Fig. 2.— Merger trees for the three most massive stars at the end 
of the simulation for the model presented in Figure Q Red shows 
the most massive, intermediate green the second most, and blue 
the third most. For each collision with the massive star, a symbol 
is plotted at its mass at the time of the collision with lines drawn 
connecting to the symbols representing the star participating in 
the merger. A triangle represents a single—single collision, a square 
a collision in a binary—single interaction, a pentagon a collision 
in a binary-binary interaction, and a circle a star that has not 
yet undergone a collision. As an example, the three intermediate 
green points at t fs 0.2 Myr represent a collision in a binary—single 
interaction of a IOOMq star with a r^O.SMq star, each having 
never undergone a collision previously. For this model there is 
clearly just one VMS, created primarily in single—single collisions. 


3 Myr. We perform simulations with binary fractions^ 
up to 0.2. The binary population is created from a clus¬ 
ter of single stars by adding secondary companions to 
randomly chosen cluster stars, with the secondary mass 
chosen uniformly in the binary mass ratio. The binary 
binding energy is distributed uniformly in the logarithm, 
truncated at high energy so that the binary members do 
not make contact at pericenter, and truncated at low en¬ 
ergy so that the orbital speed of the lightest member in 
the bina ry is larger than the local stellar velocity disper¬ 
sion Isee lFregeau et al.]l200,^ . The eccentricity is chosen 
according to a thermal distribution, truncated at large e 
so that the binary members do not make contact at peri¬ 
center. We use either N = 5 x 10® or 10® total cluster 
objects, finding no difference in the runaway results be¬ 
tween the two. We run all simulations until 3 Myr. Our 
criterion for a runaway is that its final mass is > SOOMq. 
In our runs the final masses of the VMSs are always at 
least twice this value. 

Figure ^ shows the evolution of the cluster Lagrange 
radii and the mass of the (single) runaway as a function 
of time for a cluster with tec < 3 Myr and ft = 0.05. 
The evolution is typical for models with /b < 0.1 and 
tec < 3 Myr in that: 1) the energy generated in bi¬ 
nary interactions is not sufficient to postpone core col¬ 
lapse beyond 3 Myr, and 2) there is a single runaway. 
In the model shown in Figure ^ the runaway grows to 
«28OOM0 before the cluster begins to expand in re¬ 
sponse to the energy generated in collisions with the run¬ 
away. 

We follow the membership of the collision products 

^ The binary fraction is defined such that /(, is the fraction of 
objects in the cluster that are binaries, an object being either a 
binary or a single star. 


Fig. 3.— Same as Figure 1^ for the same model but with /{, = 
0.1. A binary interaction-induced collisional runaway begins at 
t ^ 1.5 Myr and proceeds to Ril.3 X 10® M©. A direct collision- 
induced runaway begins at t ss 2.3 Myr, yielding a second runaway, 
of mass Ri2.5 X 10® Mq. The third most massive star at the end 
of the simulation has only grown to fs 4 OOM 0 . 


in binaries during the collisions and binary interactions. 
This allows us to produce merger trees to track the con¬ 
tributions to the formation of VMSs. Figure |21 shows 
the merger trees for the three most massive stars at the 
end of the simulation. There is clearly only one runaway 
for this model, and it grows primarily via direct, single¬ 
single collisions. Figure 0 shows the same as Figure |5| 
for the same model but with /& = 0.1. In this model 
there are two VMSs at the end of the simulation. A 
binary interaction-induced collisional runaway begins at 
t « 1.5 Myr and proceeds to wl.3 x 10® M©. A direct 
collision-induced runaway begins at t w 2.3 Myr, yield¬ 
ing a second runaway, of mass w 2.5 x 10® M©. The third 
most massive star at the end of the simulation has only 
grown to ~400M©. The evolution is typical for models 
with /b ^ 0.1 and tec < 3 Myr in that there are two run¬ 
aways: one via direct collisions, and one via collisions in 
binary interactions. 

Although we have not yet conducted a full param¬ 
eter space survey, we have mapped a significant slice 
in /b-trc(O) space, from /b = 0.02 to 0.2 and from 
tcc/(3Myr) = 0.16 to 28. Without exception, and in¬ 
dependent of /b, models with tec = 0.15ti.c(0) > 3 Myr 
show no runaways (7 models), while those with tec = 
0.15trc(0) < 3Myr show either one or two (10 models). 
This is in agreement with the results for clusters with 
no primordial binaries, presumably since binary fractions 
smaller than 20% are not sufficient to postpone core col¬ 
lapse beyond 3 Myr. Of the models that show runaways, 
those with /b < 0.1 always yield only single runaways 
(4 models), while those with /b > 0.1 always yield dou¬ 
ble runaways (6 models). No models ever produce more 
than two runaways. Due to the computational cost in¬ 
volved, we have not yet explored the ft, > 0.2 region of 
parameter space. 

The numerical results presented here clearly agree with 
the analytical arguments made above on the runaway 
nature of binary interaction-induced collisions. Binary 
interaction-induced runaways occur in the outskirts or 
just outside the core (typically a few percents of a parsec 
for our models). This is a region where densities, in 




























































































4 


Gurkan, et al. 


particular the binary densities, are high enough to start 
a runaway but the relaxation time is long enough that 
massive objects do not rapidly sink to the center. The 
rapid collapse near the center leads to an expansion of 
these regions, further increasing the relaxation time and 
prevents a prompt merger of the VMSs. 

Finally, we note the a pparent disagreement betwee n 
our results and those of iPortegies Zwart et alJ (1200411 . 
who performed two direct iV-body integrations of mod¬ 
els of the cluster MGG-11 with fb = 0.1 and found 
only single runaways. Our results predict that whenever 
fb ^ 0.1, and the cluster central relaxation time is short 
enough so that a runaway is expected in a corresponding 
cluster with only singl e stars, a double runaway sh ould 
result. The models of IPortegies Zwart et al.l (1200411 are 
extremely centrally concentrated, with Wq = 12, while 
our models all use Wq = 3. From the discussion above, 
it is clear that if the density profile is sufficiently steep, 
the binary interaction-induced runaway will take place 
close enough to the center so that it will seed the subse¬ 
quent direct collision-induced runaway, yielding just one 
runaway. In addition, the condition fb > 0.1 for dou- 
ble runaways is only approx imate. Since the models of 
IPortegies Zwart et al.l (1200411 are right on this boundary, 
it is possible that a small difference in the initial condi¬ 
tions has prevented a double runaway in their models. 

3. DISCUSSION 

Although the process is somewhat uncertain, it is likely 
that the VMSs formed in young clusters via collisional 
runaways will undergo core-collapse super novae and be¬ 
come IMBHs on a timescale of ~ 1 Myr (iFreitag et al.l 
l2005all . Our results show no evidence of the VMSs merg¬ 
ing prior to becoming IMBHs. After their separate for¬ 
mation, the two IMBHs will quickly exchange into a 
common binary via dynamical interactions. The IMBH- 
IMBH binary will then shrink via dynamical friction due 
to the cluster stars, to the point at which the stellar mass 
enclosed in the binary is comparable to the binary mass. 
This occurs on a timescale ~tr(w)/MiMBH, where tr is 
the local relaxation time and (m) is the local average 
stellar mass. Since (m)/MiMBH ^ 10“^, this timescale 
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is likely to be < 10 Myr. The binary will then shrink 
via dynamical encounters with cluster stars, the rate of 
which is governed by loss-cone physics. The timescale 
for the binary to shrink to the point at which it merges 
quickly via gravitational radiation energy loss is likely to 
be < 1 Gyr (jYii V Trema,inFl|2nn3t lMi]1eril2flfliTjl . 

Although IMBH-IMBH binaries do not merge in the 
LISA band—the gravitational wave frequency at merger 
is ~ 1 Hz—they do represent bright sources that take at 
least ~ 10® yr to cross the LISA band. Their inspiral 
(chirp) signals should be easily detectable by LISA out to 
a few tens of Mpc. Thus the number of detectable IMBH 
binary sources may be quite large, since most clusters are 
probably born with fb > 0.1, and any cluster with mass 
> 1O®M0 and central relaxation time < 20 Myr will lead 
to a double runaway. 

Significant core rotation (with rotational speed com¬ 
parable to the local velocity dispersion) is observed 
^in the clusters M15, Gen, 47 Tuc, and G 1 (e.g. 

rotation suggests the presence of a core angular mo¬ 
ment um source, such as an IMBH binary (iMaoelli et alJ 
[2003). Similarly, observations of a millisecond pulsar in 
the halo of NGG 6752, and two others in the core with 
high negative spin derivati ves, hint at the exi stence of an 
IMBH binary in the core (jGoIni et a,]Jl2nn3jf . Since the 
IMBH-IMBH binaries formed via collisional runaways 
will merge within 1 Gyr after formation, any angular 
momentum imparted to the cluster by the IMBH-IMBH 
binary will quickly diffuse out of the core on a core relax¬ 
ation time. An alternate mechanism must be at work in 
creating the core rotation seen in some globular clusters 
today. 
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